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Abstract — Numerical method is used in the semiconductor 
industry for the analysis and design of metal organic chemical 
vapor deposition reactors. The most important issue needed for 
prediction of epitaxial thin film growth rates and uniformity is 
the chemistry mechanism occurring in the gas phase and at the 
surface. Three numerical mechanism models are presented for 
verifying the growth rate of the gallium nitride (GaN) 
mechanism, and their results are compared with the growth 
rate results of previous literature studies. In order to reduce the 
calculation time of numerical analysis in commercial software, a 
numerical procedure is performed for simplifying the 
complicated mechanism of an epitaxial thin-film process 
through rate of production analysis. All the solution results 
could be compared in one schematic diagram, and the 
differences among these three mechanisms are extremely 
pronounced in high temperatures. The most simplified 
mechanism can be used to predict the growth rate of a 
blue-LED process before an actual process begins and it can be 
easily calculated in a standard computer. This study provides a 
useful and a predictable mechanism for GaN thin film research 
in blue-LED process. More verification work is under way by 
our own experiments. 

Index Terms — MOCVD, LED, mechanism, GaN, numerical, 
thin film. 

I. Introduction 

Metal organic chemical vapor deposition (MOCVD) is 
the key process [1~4] for manufacturing III-V compound 
semiconductor devices. III-V thin films are deposited on 
solid substrates from the gas inlet to the gas outlet. During 
this process, the reactant vapors enter the reactor from 
showerhead or nozzle and flow down vertically or 
horizontally into the reactor. Upon contacting the surface of 
the susceptor, the reactant vapors flow, radially outward 
above the susceptor and wafer then finally exhaust through 
the outlet. By using this unique technique, a high-quality 
GaN epitaxial thin film was obtained for the first time on a 
sapphire wafer [5—7] . 

Numerical modeling is now widely used in the 
semiconductor industry in the troubleshooting of epitaxial 
growth processes [8—9] . Moreover, it is commonly employed 
in key component design. In the past 20 years, only a few 
studies have investigated the GaN mechanism; therefore, the 
experimental data available for III-V nitride MOCVD growth 
chemistry are extremely limited, particularly the chemical 
kinetic rate data. Although a common reaction model is 
introduced by S. A. Safvi [10], several features of the GaN 
growth chemistry can be identified that have been established 
by experimental and numerical computational studies 


[10-14], but the GaN thin film mechanism is still unclear for 
many years. 

A complicated chemistry method derived from the ab 
initio quantum chemistry methods developed in 2005 is used 
for growth rate prediction [15]. Such detailed reaction 
mechanisms, comprising mostly elementary and ring 
reactions, have the advantage of requiring little or no 
calibration by using experimental data. However, for a 
complicated mechanism, a single computer has difficulty in 
rapidly obtaining a solution, and the calculation time is 
considerably extensive. 

For a complicate mechanism, although the accuracy is 
sufficient, the computational capacity is extremely high, the 
calculation time is always up to few days and difficult to 
calculate using a standard computer. By using the parallel 
computing may be the best way to solve this problem, but the 
fund cost is too high to achieve. This study used the 
complicated chemical mechanism from previous studies [15], 
which is simplified through the rate of production (ROP) 
analysis (CHEMKIN, 0-dimensional calculation), and the 
goal is to develop a simplified mechanism can easy 
calculated by a standard computer to predict the GaN growth 
in blue-LED process. And it can be used as the key 
equipment design for in general academia and industry. 

II. The numerical modeling details 

a. Governing equation in the 0-dimensional 

analysis 

In CHEMKIN, the Perfectly Stirred Reactor (PSR) model is 
used in this study. The contents of a PSR are assumed to be 
nearly spatially uniform due to high diffusion rates of forced 
mixing. In other words, the rate of conversion of reactants to 
products is controlled by chemical reaction rates and not by 
mixing processes. Thus we consider that the reactor is limited 
by reaction kinetics. Mass transport to the reactor walls is 
assumed to be infinitely fast. Therefore, the relative 
importance of surface reactions to gas -phase reactions is 
determined only by the surface to volume ratios of each 
material and the relative reaction rates. In this section, we 
introduce the equation that be used in this study. 

Global mass balance equation 
dm 

df. ^Hn 'dlout dl reaction (I) 

m is the mass in the reactor, rh in out is the mass flows in and 
out to the system. rh reaction is the reaction term includes 
generation and consumption work in both gas and surface 
phase. 
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Rate of production (ROP) analysis 


In CHEMKIN, ROP analysis can help user determine the 
most dominant reaction paths in some complicated reaction 
mechanisms, the general expression for this analysis is 


ROPi 


R; 


production 


i. production 




j production 


( 2 ) 


ROP, 


Ri, 


consumption 


consumption 


SjRj, 


(3) 


consumption 


ROP analysis enables the rapid identification of dominant 
reaction paths. It indicates each path by percentage on the 
x-axis (production or consumption), for example, 0.1=10%, 
0.2=20%, respectively. The y-axis shows the reaction detail 
for different temperature. Ri. production is production rate for 
the species i, £j Reproduction is the total production rate for 
species i, consumption is consumption rate for the species i, 
2j Rj, consumption is the total consumption rate for the species i. 


b. Procedure for the numerical modeling 

Fig. 1 illustrates the simplified procedure for this numerical 
modeling. The first step is 0-D computation, including an 
ROP optimized analysis for the mechanism model 
(Complete model, Reduced model, and New model) and 
general theory. 



Fig. 1 . Schematic procedure of the Complete numerical 

modeling. 


III. Results and discussion 

A. Computation of the Complete mechanism model 

A mathematical representation is developed for the Complete 
mechanism model involving ROP analysis through the 0D 
dimension software CHEMKIN. Tables 1 and 2 summarize 
the gas and surface reaction mechanisms [15]. The Arrhenius 
coefficients and activation energy are all temperature 
sensitive; the small pressure difference will not affect the 
final conclusion and computation. So, the pressure difference 
between derivative and process is ignoring. 

The gas reaction mechanism considers several aspects 
including the pyrolysis reactions for tri-methyl gallium 
(TMG), di-methyl gallium (DMG), and monomethyl gallium 
(MMG), and the adduct reactions with NH 3 and CH 3 ; 
interaction among CH 3 , H 2 , and H; and finally, the pyrolysis 
reaction for adducts. TMG and MMG decompose slowly 
because of their large activation energy barriers; however, the 
pyrolysis reaction of DMG to MMG is rapid. Therefore, the 
conversion of TMG to MMG is rapid at high temperatures, 
and the concentration of MMG is expected to be high at 
higher temperatures. Therefore, MMG and MMG:NH 3 are 
expected to be the dominant gallium-containing precursors at 


high temperatures, whereas TMG and TMG:NH 3 are 
expected to contribute to the deposition at low temperatures. 
For the surface reaction, Table 2 lists three paths and Table 3 
is the compositions for chemical compounds on the surface. 

Path 1 displays the ring reaction between MMG and NH 3 . 
First, the MMG in the gas phase is absorbed by the substrate 
to form the surface site MMG(S) which then interacts with 
NH 3 in the gas region and forms COMPMl(S). This 
procedure is repeated until COMPM5(S) is formed. 
COMPM5(S) is desorbed by a CH 3 radical to form a ring 
structure RINGM1, which is then attached to a nearby 
surface Ga(S). Finally, the hydrogen atoms H 2 are eliminated, 
resulting in GaN formation. 

Path 2 for GaN formation starts from the adsorption at the 
surface site GaNH 2 (S). This species can directly be formed 
from the gas phase reaction of the MMG and NH 3 or through 
the reaction of the adsorbed Ga(S) and adsorbed NH 2 (S). 
GaNH 2 (S) then absorbs the MMG and NH 3 in the gas phase 
to form the compound COMPMMl(S). Similar to Path 1, this 
process is repeated until the elimination of CH 4 to produce a 
ring structure RINGMl(S). Finally, the ring absorbs CH 3 
radicals from each Ga atom of the ring followed by the 
elimination of CH 4 , resulting in GaN formation. 

In Path 3, TMG adsorbs on N(S) to form TMG(s). NH 3 is 
then adsorbed by TMG(S) to form TCOMl(S). Moreover, 
TCOM1 can be formed by the adduct TMG:NH 3 in the gas 
phase, which is directly absorbed by the surface N(S). Similar 
to the processes in Paths 1 and 2, this procedure repeats until 
a complex TCOM3(S) is formed. Finally, this complex 
TCOM3(S) eliminates two molecules of CH 4 to produce 
GaN. To verify the accuracy of our CHEMKIN software, we 
incorporated the Complete mechanism model presented in 
Tables 1 and 2 into the solver. Table 4 lists the input 
parameters for the calculation of the growth rate. Fig. 2 
illustrates the predicted result compared with the results of 
previous studies. 



Temperature (K) 

Fig. 2. Result of the Complete mechanism model vs. 

literature data [11]. 

On the basis of the aforementioned mechanisms, we observe 
that in a high-temperature zone, this numerical result is 
precisely matched to the experimental data. This section 
displays that the numerical code CHEMKIN is a reliable 
software tool, and we used it to describe chemical 
mechanisms in this study. We describe the simplified 
procedure in the next section. 
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Table 1 

Gas-phase mechanism for the Complete model [15]. The unit 
of Ea is kJ/mol. 


k=AT n e Ea/RT 
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nh 3 
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CH 4 

ch 4 
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1.00E+47 

7.67E+43 

E68E+30 

2.28E+34 

E70E+04 

4.08E+31 

5.30E+05 

7.95E+24 

8.10E+05 

3.31E+03 

E20E+12 

5.00E+13 

5.00E+13 

E33E+44 

2.40E+22 

2.00E+13 

2.00E+16 


-9.18 

-9.8 

-5.07 

-8.31 

2 

-7.03 

1.56 

-5.21 

1.3 

2.51 

0 

0 

0 

-8.24 

-1 

0 

0 


76996 

34017 

84030 

3115 

19969 

3234 

20744 

2094 

17722 

9859 

12518 

10036 

10036 

77791 

0 

0 

0 


Table 2 

Surface-phase mechanism of the Complete model [15]. The 
unit of Ea is kJ/mol. 


COMPMl(S) 

COMPM2(S) 

COMPM3(S) 

COMPM4(S) 

COMPM5(S) 

RINGMl(S) 

RINGM2(S) 

COMPMMl(S) 


NH3-MMG(S) 

Ga-NH2-MMG(S) 

NH3 • Ga • NH2 ■ MMG(S) 

Ga • NH2 • Ga • NH2 ■ MMG(S) 
NH3 • Ga • NH2 • Ga • NH2 • MMG(S) 
NH2 • Ga • NH2 • Ga • NH2 • Ga( S) 
(S)NH2 • Ga ■ NH2 ■ Ga ■ NH2 ■ Ga(S) 
MMG-GaNH2(S) 


COMPMM2(S) 
COMPMM3(S) 
COMPMM4(S) 
COMPMM5(S) 
TCOMl(S) 
TCOM2(S) 
TCOM3(3) 
COMl(S) 


NH3 • MMG • GaNH2 ■ Ga(S) 
NH2 • Ga • NH2 • Ga( S) 

MMG • NH2 • Ga • NH2 • Ga( S) 
NH3 • MMG • NH2 • Ga • NH2 • Ga( S) 
NH3-TMG(S) 
NH2-DMG(S) 
(S)NH2-DMG(S) 
RINGM2(S) CH3 complex 


Table 4 

Process parameters for the 0-D model [11]. 

Volume flow rate 

values 

Standard liter per minute 

TMGa flow rate (|urnol/min) 

7.48 x 10' 4 

NH 3 flow rate (Sim) 

6 

H 2 flow rate (Sim) 

7.1 

Parameter 

value 

Initial pressure (torr) 

140 

Susceptor temperature (K) 

600 ~ 1400 

Reactor initial temperature (K) 

300 

Gas inlet temperature (K) 

300 

V/III ration 

2000 




Path 1, k 

=AT 

n.g-Ea/RT 



A 

n 

Ea 

l 

MMG + 

N(S) 


MMG(S) 



1.16E+05 

2.98 

0 

2 

MMG(S) -> 

MMG 

+ 

N(S) 



1.12E+14 

0.55 

107673 

3 

NH 3 + 

MMG(S) 

-*■ 

COMPMl(S) 



3.35E+07 

3.33 

0 

4 

COMPMl(S) -* 

nh 3 

+ 

MMG(S) 



5.70E+13 

-0.16 

8146 

5 

MMG + 

COMPMl(S) 


CH4 + 

COMPM2(S) 

1.23E+10 

3.22 

23446 

6 

NH 3 + 

COMPM2(S) 

-> 

COMPM3(S) 



3.35E+07 

3.33 

0 

7 

COMPM3(S) -> 

nh 3 

+ 

COMPM2(S) 



5.70E+13 

-0.161 

8146 

8 

MMG + 

COMPM3(S) 

-> 

CR, + 

COMPM4(S) 

1.23E+10 

3.22 

23446 

9 

NH 3 + 

COMPM4(S) 

-> 

COMPM5(S) 



3.35E+07 

3.33 

0 

10 

COMPM5(S) -> 

nh 3 

+ 

COMPM4(S) 



5.70E+13 

-0.161 

8146 

11 

COMPM51S) -> 

CH4 

+ 

RINGMl(S) 



1.23E+07 

3.22 

23446 

12 

Ga(S) + 

RINGMl(S) 

-> 

RINGM2(S) + 


N(S) 

3.35E+07 

3.33 

0 

13 

RINGM2(S) 

3H 2 

+ 

3GaN(B) + 


Ga(S) 

3.68E+09 

2.05 

59610 



Path 2, k 

=AT 

n^-Ea/RT 



A 

n 

Ea 

14 

CH3 + 

Ga(S) 

-» 

MMG(S) 



1.76E+09 

1.39 

0 

15 

MMG(S) -> 

CH3 

+ 

Ga(S) 



4.54E+13 

0.0346 

79480 

16 

NH2 + 

Ga(S) 

-> 

NH2(S) 



3.17E+08 

1.83 

0 

17 

GaNH2 + 

N(S) 

-> 

GaNH2(s) 



2.27E+06 

2.247 

0 

18 

GaNH2(S) ->■ 

GaNH2 

+ 


N(S) 



4.83E+13 

0.614 

83881 

19 

COMPMMl(S) -> 

CH4 

+ 

GaNH2(S) 



1.49E+1 1 

0.609 

25950 

20 

MMG + 

GaNH2(S) 

-> 

COMPMMl(S) 



1.16E+05 

2.98 

0 

21 

NH3 + 

COMPMMl(S) 

-> 

COMPMM2(S) 



3.35E+07 

3.33 

0 

22 

COMPMM2(S) -> 

CH4 

+ 

COMPMM3(S) 



1.49E+1 1 

0.609 

25950 

23 

MMG + 

COMPMM3(S) 


COMPMM4(S) 



1.16E+05 

2.98 

0 

24 

NH3 + 

COMPMM4(S) 

-> 

COMPMM5(S) 



3.35E+07 

3.33 

0 

25 

COMPMM5(S) -> 

CHi 

+ 

RINGMl(S) 



1.49E+1 1 

0.609 

25950 

26 

NH 2 (S) -► 

nh 2 

+ 

Ga(S) 



1.45E+14 

0.09 

59786 

27 

COMPMMl(S) -> 

MMG 

+ 

GaNH 2 (S) 



1.00E+14 

0.55 

42819 

28 

COMPMM2(S) -> 

nh 3 

+ 

COMPMMl(S) 



5.70E+13 

- 0.1 

8146 

29 

COMPMM4(S) -> 

MMG 

+ 

COMPMM3(S) 



1.00E+14 

0.55 

42819 

30 

COMPMM5(S) -> 

nh 3 

+ 

COMPMM4(S) 



5.70E+13 

- 0.1 

8146 

31 

Ga + 

N(S) 

-> 

Ga(S) 



1.00E+1 1 

1.5 

0 

32 

Ga(S) + 

NH 2 (S) 

-> 

GaNH, 

+ 

Ga(S) 

1.00E+25 

0 

0 

33 

Ga(S) -> 

Ga 

+ 


N(S) 



1.00E+13 

0 

45168 

34 

6CH; + 

RINGM2(S) 

-> 

COMl(S) 



7.55E+07 

2.31 

0 

35 

COMl(S) -> 

6 CH 3 

+ 

RINGM2(S) 



1.00E+13 

0.71 

45506 

36 

COMl(S) 

6 CH 4 

+ 

3GaN(B) 

+ 

Ga(S) 

4.00E+12 

0 

49675 



Path 3, k=AT n c 

yEa/RT 




A 

n 

Ea 

37 

TMG + 

N(S) 


TMG(S) 



1.16E+05 

2.98 

0 

38 

NH 3 + 

TMG(S) 


TCOMl(S) 



3.35E+07 

3.33 

0 

39 

TCOMl(S) ->■ 

CH 4 

+ 

TCOM2(S) 



1.49E+1 1 

0.609 

32785 

40 

Ga(S) + 

TCOM2(S) 


TCOM31S) 

+ 

N(S) 

3.35E+07 

3.33 

0 

41 

TCOM3(S) ->■ 

2 CR, 

+ 

GaN(B) 

+ 

Ga(S) 

1.49E+1 1 

0.609 

49675 

42 

TMG(S) -> 

TMG 

+ 


N(S) 



1.12E+14 

0.55 

49675 

43 

TCOMl(S) -> 

NH, 

+ 

TMG(S) 



5.70E+13 

-0.161 

11922 

44 

TMG:NH 3 + 

N(S) 

-> 

TCOMl(S) 



1.16E+05 

2.98 

0 

45 

TCOMl(S) ->■ 

TMG:NH 3 

+ 


N(S) 



1.12E+14 

0.55 

49675 

46 

TCOMl(S) ->■ 

2CH 3 

+ 

MMG(S) 

+ 

NH 3 +N(S) 

1.12E+14 

0.55 

107673 

47 

MMGNH 3 + 

N(S) 

— 

COMPMl(S) 



1.16E+05 

2.98 

0 

48 

COMPMl(S) -> 

MMG:NH 3 

+ 


N(S) 



1.12E+14 

0.55 

107673 

49 

MMG:NH 3 + 

COMPMl(S) 

-> 


CRt 

+ 

COMPM3(S) 

1.23E+10 

3.22 

23446 

50 

MMGiNH, + 

COMPM3(S) 

-> 


CH 4 

+ 

COMPM5(S) 

1.23E+10 

3.22 

23446 

51 

MMGiNH, + 

GaNH 2 (S) 


COMPMM2(S) 



1.16E+05 

2.98 

0 

52 

MMG:NH 3 + 

COMPMM3(S) 

-> 

COMPMM5(S) 



1.16E+05 

2.98 

0 

Table 3 










The chemical compositions for the chemical compounds 

on 

the surface [15] 









Compounds names 

Chemical formula 


c. Development of the Complete mechanism 

model 

For the LED industry, rapidly and accurately predicting 
'epitaxial uniformity and growth rate is critical. The accuracy 
is verified using the aforementioned 0D calculation. For the 
Complete model, the reaction number of the Complete 
mechanism is up to 69, as listed in Tables 1 and 2. Therefore, 
we must simplify the Complete model. ROP analysis is used 
for determining the contribution of each reaction to the net 
production or destruction rate of a species. It is particularly 
useful for 0D systems, where the computational expense for 
the added calculations is low, and the data can be considered 
from a large reactor set. Therefore, ROP analysis rapidly 
identifies dominant reaction paths. Fig. 3 displays the ROP 
analysis of a Complete model and the reaction path at 
different temperatures for producing a GaN thin film. For a 
blue-LED process, the critical process temperature are 
1300-1400 K. 


HOOK 

1300K 

1200K 

HOOK 

1000K 

900K 

HOOK 


600-700K 


TCOM3(S)s>2CH4+GaN(B)+Ga(S) I 
RINGM2($)=>3H2HGaN(B)+Ga(S) 
TCOM3(S)=>2CH4+GaN(B)+Ga(S) jj 
RlNGM2(S)=>3H2*3GaN(B)+Ga(S) 
TCOM3(S)=>2CH4+GaN(B)+Ga(S) 
RINGM2(S)s>3H2+3GaN(B)+Ga(S) 
TCOM3(S)=>2CH4+GaN(B)+Ga(S) 
RINGM2(S)=>3H2+3GaN(B)+Ga(S) 
TCOM3(S)=>2CH4+GaN(B)+Ga(S) 
RINGM2(S)=>3H2+3GaN(B)+Ga(S) 
TCOM3(S) s >2CH4+GaN(B)+Ga(S) 
RINGM2(S)=>3H2+3GaN(B)+Ga(S) 
TCOM3|S)=>2CH4+GaN(B)+Ga(S) 
RINGM2(S)=>3H2*3GaN(B)+Ga(S) 
TCOM3(S)=>2CH4+GaN(B)+Ga(S) 
TCOM3(S)=>2CH4+GaN(B)+Ga(S) 



Fig. 3. Rate of production analysis for the Complete model. 

For a process temperature greater than 1300 K, the most 
dominant reaction paths are reaction number 13 and 41 in 
table 1 , respectively. The percentage of the pyrolysis reaction 
for RINGM2(S) is extremely significant and can be up to 
90%. Fig. 4 depicts the schematics for Paths 1 and 2. Path 1 
displays the successive reaction between MMG and NH 3 . 
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Fig. 4. Schematic for Paths 1 and 2. 

The MMG in the gas phase is absorbed by the substrate to 
form the surface site MMG(S) and then interacts with NH 3 
until the ring structure RINGM2 is formed. The hydrogen 
atoms H 2 are then eliminated, resulting in GaN formation. 
Similar to that in Path 1, the process in Path 2 repeats until 
eliminating CH 4 ; a ring structure RINGM2(S) is then 
produced; however, the source for Path 2 is GaNH 2 . The ROP 
analysis in Fig. 4 demonstrates that at high temperatures 
(1300-1400 K), the reaction for GaNH 2 in Path 2 almost does 
not occur. 

Figs. 5, 6, and 7 depict the ROP analyses for species 
consumption of TMG, DMG, and MMG, respectively. 


1400K 

HOOK 
1200K 
HOOK 

1000K 

900K 

600 - 800K 

- 1.0 - 0.8 - 0.6 - 0.4 - 0.2 0.0 




TMG + NH, - TMGNHj 


TMG + N(S) — TMG(S) 



TMG- DMG -CHj 


HOOK 

1300K 

1200K 

HOOK 

1000K 

900K - 


600 - 800 K 



- 1.0 - 0.8 
| MMG +NH, — MMGNH, 


-0.6 -0.4 

Ring reaction 


- 0.2 0.0 
MMG + N(S) — MMG(S) 


Fig. 7. ROP analysis for MMGa in the gas phase. 

Figs. 5 and 6 display three consumption paths for TMG 
and DMG in the gas phase; however, at high temperatures, 
TMG tends to convert to DMG. A low percentage of TMG 
converts to TMG(S), and DMG tends to decompose into 
MMG in high-temperature zones. Moreover, as displayed in 
Fig. 7, MMG is consumed in the ring reaction or converts to 
MMG(S) at the surface. To sum up, we can conclude that 
nearly no CH 3 exists in the gas phase for Path 2. Thus, Path 2 
is not crucial in high-temperature zones. From [15], the 
Upper way reaction (ring structure) is the main way to 
produce the GaN, every pathway is first connected by three 
Ga sources and three N sources, then it connected by another 
active point, finally release the CH 4 to produce the GaN, but 
this is too complicate to calculate. So, in this section, this 
study assumed that a single direction way to produce GaN, 
when the Ga source adsorbed by the surface, the NH 3 in the 
gas phase will physically adsorbed on the surface, when the 
energy is sufficient, it will release CH 4 , H 2 to form GaN. 

Hence, we delete Path 2 and the ring structure for Path 1 
for simplifying the Complete model. Fig. 8 displays a 
schematic of the new Path 1 (reduced surface mechanisms). 


Change the Path 1 


Reduce Model 


TMG 


MMG 


: r. 
c 


NH, 





COMPM 1 


ch 4 






a 

1 

i 

MMG:NHj 



L-»| RINGMI 



Fig. 5. ROP analysis for TMGa in the gas phase. 


Fig. 8. Schematic diagram for the new Path 1. 


1400K 

HOOK 

1200K 

HOOK 

1000K - 



-08 -0.6 
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Fig. 6. ROP analysis for DMGa in the gas phase. 


Table 5 lists the new form of surface mechanisms in 
Reduced model (gas mechanism remains the same from 
complicate mechanism) for Path 1 . 


Table 5 


Surface mechanisms (Path 1+Path 3) are considered for the 


Reduced model. 


1 MMG + N(S) MMG(S) 

2 MMG(S) -> MMG + N(S) 

3 NH 3 + MMG(S) COMPMl(S) 

4 COMPMl(S) -> NH 3 + MMG(S) 

5 COMPMl(S) -> CH 4 + COMPM2(S) 

6 COMPM2(S) -* H 2 + GaN(B) + N(S) 

7 TMG + N(S) -> TMG(S) 

8 TMG(S) -> TMG + N(S) 

9 NH 3 + TMG(S) -*■ TCOMl(S) 

10 TCOMl(S) -* NH 3 + TMG(S) 

11 TCOMl(S) - CH 4 + TCOM2(S) 

12 TCOM2(S) - 2CH 4 + GaN(B) + N(S) 


A 

n 

Ea 

1.16E+05 
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0 

1.12E+14 

0.55 

107673 

3.35E+07 

3.33 

0 

5.70E+13 
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Fig. 9 provides the final predicted result (include path 1 
and path 3) of this new form of mechanism for the Reduced 
model; the result strongly agrees with that of the Complete 
model. 


3.0 


25 


£ 2.0 

co 

c 

0 

2 i( . 

1 15 

£ 

(T 10 
si 

S 

2 o 5 


0.0 - • •- 


600 


* 


■ 


• Complete model 

• Reduced model 


800 


1000 
T (K) 


1200 


1400 


Fig. 9. Complete model vs. the Reduced model. 


d. Mechanism optimization 

We investigated whether the mechanism (Reduced model) 
can be simplified further. As illustrated in Fig. 6, in a 
low-temperature zone, TMG tends to form an adduct with 
NH 3 . However, in a high-temperature zone, all the TMG is 
rapidly consumed by DMG. Moreover, as displayed in Fig. 7, 
similar to TMG, in a high -temperature zone, all the DMG is 
consumed by MMG; however, DMG tends to form an adduct 
with NH 3 in a low-temperature area. However, as illustrated 
in Fig 8, at a high temperature, MMG participates in the ring 
reaction at the wafer surface; therefore, for a blue-LED 
process, we can perform the next simplified procedure. 


Fig. 10 provides the final predicted results of the Complete 
model, Reduced model, and New model, and the 
experimental data of their comparisons. The difference 
between the New model and the Complete model is low. 



Fig. 10. Comparison among all mechanism models. 


The pressure effect of growth rate for New model 

Fig. 1 1 provides the pressure effect of GaN thin film growth 
rate for New model in different temperature, the pressure 
effect is significant only on 1400K. When the low 
temperature reaction, the GaN reaction mechanism is 
dominant by surface kinetic, so the pressure is not a key 
factor in low temperature. In 1400K, the growth rate is 
change when the reaction pressure lower than 200torr, this is 
contributed by the reduced for residence time. 



Fig. 1 1 . The pressure effect for the New model 


In a high-temperature zone, we assume that all TMG is 
consumed to immediately form MMG, and the adduct 
assumption is negligible. Therefore, for a gas reaction, the 
mechanism can be reduced to only one chemical reaction, 
and the reaction for adduct is negligible; Table 6 presents the 
results (New model). 


The III/V ratio effect of growth rate for New model 

Fig. 12, 13 provides the effect of NH 3 and TMGa, we can 
clearly see that the TMGa is the dominant species in the full 
New mechanism model, when the gas flow rate of the NH 3 
increase, the residence time is reduced, but this effect is lower 
than surface reaction by NH 3 . 


Table 6 


New mechanisms (New model) for GaN deposition. 




Gas Reaction, k=AT n e Ea/RT 

A n 

Ea 

1 

TMG 

MMG + C 2 H 6 

1.00E+47 1.00E+47 

-9.18 

Surface Reaction 


2 TMG + N(S) 

3 TMG(S) TMG 

4 NH 3 + TMG(S) 

5 TCOMl(S) -» NH 3 

6 TCOMl(S) 3CH 4 

7 MMG + N(S) 

8 MMG(S) MMG 

9 NH 3 + MMG(S) 

10 COMPMl(S) -> NH 3 

11 COMPMl(S) -> H 2 



TMG(S) 


1.16E+05 

2.98 
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Fig. 12. The effect of V/III ratio (The flow rate of species V is 

constant) 
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IV. Conclusions 

A simplified process that used CHEMKIN for modeling the 
GaN growth is demonstrated. Several reaction pathways are 
calculated using ROP analysis. The activation barriers and 
vibrational frequencies are employed for calculating rate 
constants. The predicted growth rates agreed well with the 
experimental data at each temperature point. The most 
simplified mechanism can be used to predict the growth rate 
of a blue-LED process before an actual process begins; this is 
extremely useful for LED companies. All of these 
observations demonstrate that the proposed mechanisms are 
fairly universal and crucial pathways for the epitaxial growth 
of GaN. This study provides a useful and a predictable 
mechanism for GaN thin film research in blue-LED process. 
More verification work is under way by our own 
experiments. 
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